Effect of wavelength dependence of nonlinearity, gain, and dispersion in photonic 

crystal fiber amplifiers 
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Photonic crystal fibers are used in fiber amplifiers and lasers because of the flexibility in the design 
of mode area and dispersion. However, these quantities depend strongly on the wavelength. The 
wavelength dependence of gain, nonlinearity and dispersion are investigated here by including the 
wavelength dependence explicitly in the nonlinear Schrodinger equation for photonic crystal fibers 
with varying periods and hole sizes. The effect of the wavelength dependence of each parameter is 
studied separately as well as combined. The wavelength dependence of the parameters is shown to 
create asymmetry to the spectrum and chirp, but to have a moderating effect on pulse broadening. 
The effect of including the wavelength dependence of nonlinearity in the simulations is demonstrated 
to be the most significant compared that of dispersion or gain. 

PACS numbers: 



INTRODUCTION 

Photonic crystal fibers are a new class of optical fibers 
that have a periodic cladding Q, 0] ■ Light can be con- 
fined to the core either by the band gap effect or by total 
internal reflection by average refractive index difference. 
Photonic crystal fibers have many intriguing character- 
istics, for example, they can be endlessly single mode, 
have extremely small or large mode areas and still re- 
main single mode. Also the dispersion properties of pho- 
tonic crystal fibers are very different from standard op- 
tical fiber. The research of photonic crystal fiber lasers 
Sl[iJi|[iJ^^ and amplifiers 

lllllilialillililijhas been intense in the past 
few years. The aim has been two-fold. First, the possibil- 
ity of obtaining a small mode area has been exploited by 
designing high-gain amplifiers/lasers where the overlap 
between the mode distribution and doped area is maxi- 
mized. On the other hand, the possibility of obtaining a 
large mode area has been utilized for realizing high-power 
amplifiers /lasers with low nonlinearity. 

Small-mode area photonic crystal fibers generally 
have large dispersion and nonlinearity that also depend 
strongly on the wavelength 25] . In this paper, we study 
the propagation of 200 fs pulses in high-gain small-mode 
area photonic crystal fiber amplifiers. When short pulses 
are considered, the spectrum is wide and thus the wave- 
length dependence of the different parameters has a pro- 
found effect on the pulse propagation. 

We compare the temporal and spectral widths, time- 
bandwidth products, and chirps of the pulses and am- 
plification properties for different fiber geometries. We 
find out that the wavelength dependence of the parame- 
ters has a substantial effect on the pulse properties after 
it has propagated a short distance in the fiber ampli- 
fier. The dispersion caused by the wavelength depen- 
dence of the nonlinearity counteracts the dispersion of 



the fiber and thus the pulses do not broaden as much 
as it is expected by considering constant parameter val- 
ues for dispersion and nonlinearity. On the other hand, 
the spectral broadening and chirping become asymmet- 
ric when the wavelength dependence of the parameters is 
taken into account. Also, the wavelength dependence of 
the dispersion parameters and gain is seen to influence 
the pulse quality less, compared to the important effect of 
wavelength dependence of nonlinearity, indicating that in 
some cases they could be approximated by constant val- 
ues when simulating pulse propagation in photonic crys- 
tal fibers. 



NUMERICAL METHODS 

The pulse propagation is studied with the optical non- 
linear Schrodinger equation 
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which is simulated by the split-step Fourier method [26j . 
The parameters for gain, nonlinearity and dispersion 
are taken to be either constants g, 7 and j3 m or wave- 
length dependent g(X), 7(A) and (3 m (A) , respectively. 
The slowly varying envelope of the pulse is taken to be 
Gaussian 



A(z = 0,T) = ^P~ e- T2 ^ 2T o\ 



(2) 



where Pq is the peak power and Tq is the pulse length. 

The values of the dispersion and nonlinear parameters 
are calculated as explained in Ref. [22]] • The dispersion 
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parameters f3 m are defined as 
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The mode-propagation constants (3 as a function of the 
frequency u> were calculated with the full-vectorial plane 
wave method (the MIT Photonic bands software) [28j |. 
The nonlinear parameter 
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is inversely proportional to the effective area, which is 
calculated from the intensity distribution of the eigen- 
mode 
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and ni is the nonlinear-index coefficient n-i 
10- 20 m 2 /W. 

We consider an Erbium-doped fiber amplifier. The 
gain as a function of wavelength is as in Fig. 1 of Ref. [13 
and it is approximated to be the same for all geome- 
tries since it has the least effect on the pulse propagation 
characteristics such as pulse broadening. Also, accord- 
ing to Ref. [3 the gain dependence on the wavelength is 
strongly influenced by the emission and absorption cross 
section rather than by the photonic crystal fiber geome- 
try. 



THE STUDIED FIBER GEOMETRIES AND 
PULSE PROPERTIES 

We investigate photonic crystal fibers with a triangular 
lattice of air holes in the cladding. The core is formed 
by a missing hole. We consider geometries with the hole 
diameter to period ratios d/P = 0.2, 0.3, 0.4, 0.5, 0.6 
and periods P = 2, 3, 4, 5, 6 fim. The geometries are 
those of interest in the research on high-gain efficiency 
photonic crystal fiber amplifiers [TBL l2ll 1221 ] . The 
dispersion and nonlinear parameters [see Eqs. © and 
(0}] are calculated as a function of the wavelength for all 
the geometries. They are shown in Fig. ^ f° r t he fiber 
geometry d/P = 0.3. 

The magnitude of the period affects the wavelength 
dependence of the (3 m . Regardless of d/P, for P = 2 
fxm, the 0m are strongly dependent on the wavelength 
while for P — 6 /im the wavelength dependence is not as 
prominent. The functional form of the m as a function 
of the wavelength is very different for the different fiber 
geometries. For example, for P = 2 /im, the value of 
fa in the center of the considered wavelength range is 
positive when d/P=0.2, 0.3, and 0.4, but negative when 
d/P=0.5 and 0.6, which can affect the pulse propagation 
considerably. The nonlinear coefficient 7(A), however, 



has a similar form as a function of wavelength for all 
fiber geometries. The steepness of 7(A) increases when 
P decreases or d/P increases. 

The studied pulse has the length 200 fs, wavelength 
1.55 /urn, and peak power P =0.01 W. The propaga- 
tion distance is 10 cm. The time and frequency axis are 
divided into 4096 steps in the split-step Fourier-method. 

To demonstrate the effect of the wavelength depen- 
dence of g(X), 7(A), and /3 m (A) in Eq. @J, the simula- 
tions for all the different geometries are performed twice: 
with wavelength dependent gain, nonlinearity and dis- 
persion parameters and, for comparison, with all these 
parameter values constant. The constant parameter val- 
ues are determined at A = 1.55 /im. The simulations for 
some of the geometries (d/P = 0.4 and d/P = 0.5) are 
repeated keeping the nonlinear parameter and gain wave- 
length dependent, but approximating the the dispersion 
parameters with constant values. Also, to investigate the 
effect of the wavelength dependence of gain, the simula- 
tions are repeated with all other parameters constant but 
including the wavelength dependence of gain. 



COMPARISON BETWEEN SIMULATIONS WITH 
CONSTANT AND WAVELENGTH DEPENDENT 
GAIN, NONLINEARITY, AND DISPERSION 

We compare the simulations made with wavelength de- 
pendent and constant gain, nonlinear and dispersion pa- 
rameters. After propagating 10 cm in the fiber amplifier, 
which is a short distance compared to actual amplifier 
lengths, there are clear differences for the two sets simu- 
lations. The pulse shape, spectrum, phase, and chirp are 
shown in Fig. for one fiber geometry. The definitions 
of the phase of the pulse <f> and frequency chirp 8ui are 



A{z,T) = \A{z,T)\eW 
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The chirp is a measure how much the instantaneous fre- 
quency changes across the pulse from the central fre- 
quency [2(|. From Fig. [21 one can see that the self- 
steepening of the pulse is larger for the simulation with 
constant g, 7, and f3 m than for wavelength dependent 
g(X), 7(A), and /3 m (A). The spectrum of the pulse is 
more asymmetric when the wavelength dependence of 
the parameters is taken into account. Also, the chirp 
of the simulation with constant parameters is symmetric 
whereas the chirp of the one with the wavelength depen- 
dent parameters is larger (in absolute magnitude) on the 
leading edge of the pulse than on the trailing edge. 

To characterize the pulses after the propagation in the 
amplifier, we calculate the temporal and spectral widths, 
time-bandwidth products, and chirps of the pulses as well 
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FIG. 1: Dispersion and nonlinear parameters as a function of 
wavelength for a photonic crystal fiber with d/P = 0.3 and 
periods P = 2 /im (blue), P = 3 //m (red), P = 4 /xm (green), 
P = 5 /im (magenta), and P = 6 ^m (black). 




FIG. 2: (a) Pulse shape, (b) spectrum, (c) phase, and (d) 
chirp after it has propagated 10 cm in the fiber amplifier. 
The fiber geometry is d/P = 0.3 and P = 3 fim. Blue curves 
represent the simulation with wavelength dependent parame- 
ters [g(X), 7(A), /3 m (A)] and red dashed curves with constant 
parameters (g, 7, fim). The propagation distance is short 
compared to the total lengths of the fiber amplifiers (which 
can be several meters). 
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as their amplification in the fiber. The temporal width is 
described by the full-width at half-maximum (FWHM) 
and denoted by AT FWHM . We also calculated the root- 
mean square widths of the pulses, but they are not shown 
here since there was no additional information compared 
to the FWHM widths. The spectral width is calculated 
as the full- width at half-maximum Acofwhm- The pulse 
widths, the time-bandwidth products 



TBP 



AT, 



FWHM 



Alj 



FWHM 



2tt 



(8) 



and the chirps are shown in Fig. |3| as a function of the 
period for the different fiber geometries. The chirp is 
indicated by the lowest value of the chirp [compare to 
Pig. Eld)]. 

The results for constant and wavelength dependent g, 
7, and f3 m are different. Pulses broaden less when the 
wavelength dependence of the parameters is taken into 
account, for all fiber geometries except d/P — 0.2. This 
indicates that the wavelength dependence of the nonlin- 
earity counteracts the temporal broadening induced by 
the dispersion. As was discussed in Section „ the form of 
the nonlinearity as a function of the wavelength is simi- 
lar for all the fiber geometries, only the steepness of 7(A) 
decreases for decreasing d/P, whereas the form of the 
dispersion as a function of the wavelength is completely 
different for all fiber geometries. From Fig. E2a) one can 
see that the temporal broadening of the pulses is reduced 
for all fiber geometries, only for d/P — 0.2 the effect has 
vanished, due to the slight steepness of the nonlinearity 
curve. 

The magnitude of the chirping is not strongly affected 
when the wavelength dependence of the parameters is 
included, but it shows a clear asymmetry that can be 
seen in Figs.EIc) and^d). The chirp magnitudes for the 
geometries with P = 2 /Ltm, d/P = 0.2,0.3,0.4,0.5,0.6 
and for P = 3 fj,m, d/P = 0.6, are not shown in Fig.OJd) 
since the chirps are linear across the total length of the 
pulse and thus exhibit no relevant minima. 

It can be seen from Fig.[3]that the pulses broaden more 
when d/P increases or P decreases, because the nonlin- 
ear and dispersion parameters increase with increasing 
d/P and/or decreasing P. For fiber geometries that have 
large periods, the d/P has less effect on the pulse broad- 
ening, because the dispersion and nonlinearity of large 
d/P geometries do not depend strongly on P. The ge- 
ometries with small P and d/P are not suitable for short 
pulse amplification since the broadening of the pulses is 
so strong. Also, the chirping lowers the quality of the 
pulses and effects dispersion compensation schemes. 

The amplification in the 10 cm of the fiber amplifier 
was calculated as 
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FIG. 3: (a) Temporal full-width at half-maximum, (b) spec- 
tral full- width at half- maximum, (c) time-bandwidth product, 
and (d) chirp of the pulses after propagating 10 cm in the 
fiber amplifier for the different fiber geometries: d/P — 0.2 
(magenta curves), d/P = 0.3 (blue curves), d/P = 0.4 (red 
curves), d/P = 0.5 (green curves), and d/P — 0.6 (black 
curves) as a function of the period P. Solid curves represent 
the simulations for which the wavelength dependence of g, 
7, and p m is taken into account. Dashed curves represent 
the simulations where these parameters are constant. The 
pulse widths are scaled by the width of the initial Gaussian 
pulse. The chirp is the lowest value of the chirp [compare to 
Fig.Hd)]. 



FIG. 4: The amplification of the pulses after propagating 10 
cm in the fiber amplifier for the different fiber geometries: 
d/P — 0.2 (magenta curve), d/P — 0.3 (blue curve), d/P = 
0.4 (red curve), d/P = 0.5 (green curve), and d/P — 0.6 
(black curve) as a function of the period P. 

and it is shown in Fig. ^ Here the wavelength depen- 
dence of the parameters is taken into account. The fibers 
with small P and/or large d/P are amplified less since the 
spectrum of the pulses broadens heavily and not all fre- 
quency components get amplified since the Erbium gain 
spectrum is so limited. This indicates that the effects of 
pulse spectrum broadening have to be taken into account 
when high-gain fiber amplifiers are designed. 

The results of two fiber geometries {d/P = 0.3, P = 2 
/j,m and d/P = 0.4, P — 2 fim) are not shown in Figs. El 
andEl These fiber geometries have very large third order 
dispersion which results in a complicated pulse shape (see 
Fig. EJ). The behavior of these pulses is completely differ- 
ent for the two simulations with constant and wavelength 
dependent parameters (compare the blue and red dashed 
curves in Fig. EJ). Generally, when P is small small, the 
parameters (3 and 7 depend strongly on the wavelength, 
and they cannot be approximated by constant values. 

COMPARISON OF THE IMPACT OF THE 
WAVELENGTH DEPENDENCE OF THE 
DIFFERENT FIBER PARAMETERS 

In order to compare the effect of including the wave- 
length dependence of the different fiber parameters on 
the pulse characteristics, we calculated the pulse propa- 
gation for the fiber geometries d/P = 0.4 and d/P = 0.5 
with two new sets of simulations. Firstly, we approxi- 
mated the dispersion parameters j3 m with constant val- 
ues, and secondly, we approximated all other parameters, 
except gain, with constant values. The time-bandwidth 
products are shown in Fig. El From Fig. Et a ) one can 
see that the results for the simulations where all the pa- 
rameters were wavelength dependent (blue squares) and 
where dispersion parameters (3 m were approximated by 
constant values (red stars) are very close, compared to 
the simulations with all parameters constant (black cir- 
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FIG. 5: (a) Pulse shape and (b) spectrum after propagating 9 
cm in the fiber amplifier with geometry d/P = 0.4 and P = 2 
fim. Blue and red dashed curves represent the simulation 
with wavelength dependent and constant parameters g, 7, 
Pm, respectively. 



cles). This indicates that the wavelength dependence of 
the dispersion coefficients /3 m does not have as signifi- 
cant effect to the pulse as the wavelength dependence of 
the nonlinearity coefficient. For the fiber geometries in 
Fig. HJb) the wavelength dependence of the dispersion 
has a larger impact on the pulse characteristics, since 
the pulse spectrum for the geometries with larger d/P is 
wider. 

To determine the effect of the wavelength dependence 
of gain, we calculated the pulse propagation when only 
the wavelength dependence of gain was taken into ac- 
count but all other parameters (/? m and 7) were constant. 
The time-bandwidth products are shown in Fig. El One 
can see in Fig. EI a ) that the simulation where the wave- 
length dependence of gain is taken into account (green 
crosses) does not change the pulse properties in a sig- 
nificant amount compared to the simulation where all 
the parameters were constant (black circles). Again in 
Fig-Elk) the effect of the wavelength dependence of gain 
is more important due to the extensive pulse broaden- 
ing. This implies that the wavelength dependence of gain 
does not have a significant effect to the pulse propagation 
except in the case where the spectral broadening of the 
pulse is excessive. For the fiber geometry with d/P = 0.5 
and P = 2 /im, it is very important to include the wave- 
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FIG. 6: Time-bandwidth products for geometries d/P = 0.4 
(a) and d/P = 0.5 (b) as a function of the period for four 
different sets of simulations. Blue squares denote simula- 
tions with all parameters wavelength dependent and red stars 
denote simulations where dispersion constants were approxi- 
mated with constant values. Black circles denote simulations 
with all parameters constant and green crosses denote simu- 
lations where only gain is wavelength dependent. The TBP 
values are scaled to the results from the simulations with all 
parameters wavelength dependent for the corresponding fiber 
geometry. 



length dependence of gain in the simulation since the 
pulse broadens out of the gain spectrum of Erbium. 



CONCLUSIONS 

We have studied pulse propagation in high-gain effi- 
ciency photonic crystal fiber amplifiers with varying pe- 
riods and hole sizes. We took into account the wavelength 
dependence of the fiber parameters for dispersion, non- 
linearity, and gain. The wavelength dependence of the 
fiber parameters has a significant effect on the temporal 
and spectral width of the pulse. The pulses were shown 
to broaden less for most fiber geometries when the wave- 
length dependence of the parameters was taken into ac- 
count, indicating that the wavelength dependence of the 
nonlinearity counteracts the dispersion of the fibers. The 
spectral width and chirp showed asymmetry after propa- 
gating a short distance in the fiber amplifier. This could 
affect for example pulse compression or dispersion com- 
pensation schemes. Although the changes in the pulse 



properties shown here were qualitatively small, they are 
important since the propagation distance was short com- 
pared to actual amplifier lengths. The wavelength de- 
pendence of dispersion and gain was shown not to have 
as profound effect on the pulse quality as the wavelength 
dependence of nonlinearity However, the wavelength de- 
pendence of all fiber parameters have to be included in 
the simulations when the pulse spectrum broadens heav- 

iiy- 
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